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Recent theoretical predictions and experimental measurements have demonstrated that equilib- 
rium free energy differences can be obtained from exponential averages of nonequilibrium work 
values. These results are similar in structure, but not equivalent, to predictions derived nearly three 
decades ago by Bochkov and Kuzovlev, which are also formulated in terms of exponential averages 
but do not involve free energy differences. In the present paper the relationship between these two 
sets of results is elucidated, then illustrated with an undergraduate-level solvable model. The anal- 
ysis also serves to clarify the physical interpretation of different definitions of work that have been 
used in the context of thermodynamic systems driven away from equilibrium. 



I. INTRODUCTION 

In recent years there has been considerable interest in the nonequilibrium statistical mechanics of small systems [l[ . 
Among the results that have been derived and tested experimentally, the nonequilibrium work theorem 0, Q , 

( e -P w ) = e-P AF , (1) 

relates fluctuations in the work W performed during a thermodynamic process in which a system is driven away from 
equilibrium, to a free energy difference AF between two equilibrium states of the system. Here, (3 specifies an inverse 
temperature, and the angular brackets denote an average over an ensemble of realizations (repetitions) of the process 
in question Q . Eq. [T] and closely related results 0, @, 0] , along with experimental confirmations d, H, [l(J H3 , have 
revealed that equilibrium free energy differences can be determined from distributions of nonequilibrium work values. 

The recent progress in this area has drawn attention to a set of earlier papers by Bochkov and Kuzovlev [l2|, fl3l . 
Hi], [HI, in which the authors had obtained - as one consequence of a more general analysis - the following result: 

(e-P w °) = 1. (2) 

The angular brackets and inverse temperature fj appearing here have the same meaning as in Eq. [IJ and Wo is 
identified as the work performed on the system. 

Although Eqs.[T]and[2]are evidently similar in structure, they are not identical; most notably, AF does not appear, 
either explicitly or implicitly 1 , in Eq. [51 The precise relationship between these two results has not been clarified 
in the literature, nor is it immediately obvious from a quick comparison of the original derivations. The aim of the 
present paper is to fill this gap, first by deriving the two equalities within a single, Hamiltonian framework, and then 
by illustrating them both using the simple model of a perturbed harmonic oscillator. The conclusions that will emerge 
from this analysis are summarized by the following three points. 

• Eqs.[T]and[2]apply to the same physical situation: a system, initially described by an unperturbed Hamiltonian 
Hq, is driven away from equilibrium by the application of a time-dependent perturbation. In principle, a single 
set of experiments could be used to test both predictions. 

• While both W and Wo are identified as work (in Refs.JH, 0] and [12, 0, 0, EH , respectively), the two quantities 
generally differ; see Eq. [751 below. The difference between them amounts to a matter of convention, related to 
whether we choose to view the perturbation as an external disturbance, or else as a time-dependent contribution 
to the internal energy of the system. 

• For the special case of cyclic processes, in which the perturbation is turned on and then off, Eqs. Q] and [2] are 
equivalent. 



Rewriting Eq.[T] in terms of dissipated work Q], = W — A_F, we obtain (exp(— /3W^)) = 1, which bears an even stronger resemblance 
to Eq.[2] However, the quantity Wq appearing in Eq. [2] is not equivalent to W^, as apparent from the definitions provided in Section ITT1 
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This paper is organized as follows. Section |TT] establishes the Hamiltonian framework and the notation that will be 
used throughout the paper. In Section Hill we derive Eqs. [1] and [2] within this framework. Section ITVl describes an 
exactly solvable model - a harmonic oscillator driven by a time-dependent external force - that illustrates the validity 
of these predictions and provides intuition regarding the two definitions of work, W and Wq. Finally, Section [Vl 
presents an alternative derivation of Eqs. Q] and [2] by way of a stronger set of results (Eq. [57)) . The paper concludes 
with a brief discussion. 



II. SETUP 



To carry out a direct comparison between Eqs. [T] and [21 we will use the setup considered in Ref. [14| . Consider 
a classical mechanical system with D degrees of freedom, described by coordinates q = (qi, ■ ■ ■ ,qo) and momenta 
p = (pi, • • ■ ,Pd), and let z — (q, p) denote a point in the phase space of this system. Consider also a number of 
external forces X\,X%, ■ ■ ■ , which are under our direct control. We act on the system by manipulating these forces. 
The Hamiltonian that describes this system takes the form 

H{z;X)=H (z)-J2 x kQk{z) (3) 

fc 

(see Eq. 2.2 of Ref. [3), where Qi(z), Q2(z), ■ ■ ■ denote the variables conjugate to the external forces: 

dH 



dx k 



(4) 



H is a function on phase space, parametrized by the forces X = (Xi, X-2, ■ ■ • ). We will refer to Hg as the bare, or 
unperturbed, Hamiltonian, and to H as the full Hamiltonian. 

If this system is brought into weak contact with a thermal reservoir at temperature T, with the external forces held 
fixed, then it will relax to an equilibrium state described by the Boltzmann-Gibbs distribution 

P cq (z;X)= * exp X)], (5) 

where (3 = (fceT) -1 . The corresponding classical partition function and free energy are: 

Z(X) = J dz exp[-/3H(z;X)] , F(X) = -/T 1 hiZ(X). (6) 

Now imagine that we subject this system to a thermodynamic process, defined by the following sequence of steps. 
Prior to time t = 0, the system is prepared in equilibrium, in the absence of external forces, i.e. at 

X = (0,0,---). (7) 

The reservoir is then removed. Subsequently, from t — to a later time t = t, the external forces are turned on 
according to some arbitrary but pre-determined schedule, or protocol, X t . The microscopic evolution of the system 
during this interval of time is described by a trajectory z t evolving under Hamilton's equations, 

dq dH dp dH 

~dt = ~dp~ ' ~dt = ~ 5q"' ( ' 

where H = H(z;X t ). The protocol X t effectively traces out a curve in "force space", from the origin (Eq. [7]) to 
some final point X T . Let AF denote the free energy difference between two equilibrium states - both at the same 
temperature T - associated with the initial and final forces: 

AF = F(X T ) F(X Q ) = -/T 1 In . (9) 

By repeatedly subjecting the system to this process - always first preparing the system in equilibrium, and always 
using the same protocol X t - we generate a number of statistically independent realizations of the process, each 
characterized by a Hamiltonian trajectory z t describing the microscopic response of the system to the externally 
imposed perturbation. Angular brackets (• ■ • } will specify an ensemble average over such realizations. 
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For a given realization, let us now define W and Wo appearing in Eqs. [T] and [2] 

W = f dtJ2x k (t)Q k (z t ) (10a) 
Jo k 

W = - [ dtJ2x k (t)Q k {zt), (10b) 
Jo k 

where the dots denote time derivatives, e.g. Q k = (d/dt)Q k (zt). These two definitions are not equivalent: in general, 
W ^ Wo- 

To gain some physical insight into these quantities, we rewrite them as follows: 

Wo = JdtX-Q = Jx-dQ (11a) 

W = -JdtX-Q = JdX~, (lib) 

where Q — (Q±, Q2, • • • ) is the vector of variables conjugate to the forces X = (X±, X 2 , ■ ■ ■ ) (sec Eq.S]). The expression 
for Wo is the familiar integral of force versus displacement found in introductory textbooks on mechanics [16j |. and 
corresponds to the definition of work used by Bochkov and Kuzovlev (Eq. 2.9 of Ref. [H|]). By contrast, expressions 
equivalent to Eq. Illbl are often used to define work in discussions of the microscopic foundations of macroscopic 
thermodynamics jl7l. Il8l. Il9j: this is the definition that is used in the context of nonequilibrium work theorems (e.g. 
Eq. 3 of Ref. 0]). While it might seem unusual that two different quantities, Wq and W , can both be interpreted as 
the work performed on a system, this ambiguity simply reflects the freedom we have to define what we mean by the 
internal energy of the system of interest. We discuss this point in some detail in the following two paragraphs. 

What is the internal energy of the system when its microstate is z — (q, p) , and the external forces are set at values 
X = (Xi, X2, • • • )? Eq. [3] suggests two natural ways to answer this question, (i) We can take the internal energy to 
be given by the value of the bare Hamiltonian, Ho(z). From this perspective the system is imagined as a particle in 
a fixed energy landscape, Ho] we affect the particle's energy by varying the forces so as to move it from one region of 
phase space to another, but the forces X do not themselves appear in the definition of its energy, (ii) Alternatively 
we can define the internal energy to be given by the value of the full Hamiltonian, H — Hq — X ■ Q. This point of view 
is captured by imagining an energy landscape that is not fixed, but changes with time as we manipulate the forces 
X. Let us refer to these two alternatives as the (i) exclusive and the (ii) inclusive frameworks, according to whether 
the term —X ■ Q is treated as a component of the internal energy of the system. 

Now we use the Hamiltonian identity 

dH ch_dH_ dH__dH_ dH _ 
dz dt dq dp dp dq 

(see Eq. [5) to obtain 

d TTI ^ \ dH dz dH dX ■ dH ■ „ 

dt H ^ x ^ = ^-Tt + ^-^ = x -dx = - x - Q > (13) 

and therefore 

j t H (z t ) = j t (H + X-Q)=X-Q. (14) 

Comparing with Eq. QTJ we see that Wo and W are equal to the net changes in the values of Ho and H, respectively, 
during the interval of perturbation: 

Wo = J dt ^ = H (z T ) - H (z ) (15a) 

W = J dt—=H(z r ]X r )-H(z ;X ). (15b) 

Since the system is thermally isolated (i.e. not in contact with a heat reservoir) from t = to t = t, it is natural 
to identify the work performed on it with the net change in its internal energy. With this in mind, Eq. [T5] provides 
a simple interpretation of the difference between W and Wo- If we adopt the exclusive point of view and take the 
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internal energy to be the value of the bare Hamiltonian Hq, then Wq is the work performed on the system, by the 
application of external forces that affect its motion in a fixed energy landscape. If we instead choose the inclusive 
framework, using the full Hamiltonian H = Hq — X ■ Q to define the internal energy of the system, then W is the 
appropriate definition of work. The distinction between these two frameworks is illustrated with a specific example 
in Section HVl 

From Eq. [15] we obtain an explicit expression for the difference between W and Wq : 

W -W = X T -Q(z r ) = J2 x k(r)Qk(zr), (16) 

k 

since X = (0, 0, • • • ). 



III. DERIVATIONS 



Let us now compute the averages of e~^ w ° and e~@ w , over an ensemble of realizations of the thermodynamic 
process described above. Since the system evolves under deterministic (Hamiltonian) equations of motion from t = 
to t = r, a given realization is uniquely determined by the initial conditions zq. We can therefore express as 
an integral over an equilibrium distribution of initial conditions: 



(e-«) = / dz P c(i (z ;X )e-^ w ^ : 



(17) 



where Wq(zo) denotes the value of Wo for the trajectory launched from the microstate zq. The first factor in the 
integrand is 



P e *(z ;X ) = 



Z(X ) 



o -pH (z ) 



(18) 



[note that H(z Q ;X ) = H {zq), by Eq. [7]- Using Eq. I15a| we have 

Wo(z ) = H a (z T {zo)) - Ho(zo), (19) 

where z T (zo) indicates the final microstate of this trajectory, expressed as an explicit function of the initial microstate. 
Upon substituting these expressions into Eq. 1171 a cancellation of terms occurs in the exponents, and we get 



(e-P Wo ) = 



Z(X ) 



dz e-P Ho(z ^o)). 



(20) 



Since there is a one-to-one correspondence between the initial and final conditions of a given trajectory, we can change 
the variables of integration from zq to z T : 



(e-P w °) 



Z(X ) 



dz T 



dz T 



dz 



-0H o (z T ) 



(21) 



We have inserted the determinant of the Jacobian matrix associated with this change of variables. By Liouville' 
theorem, this factor is identically unity, \dz T /dzo\ — 1, which finally gives us 



Z(X ) 



dzre- 01 * ^ = 1, 



(22) 



byEq.l 

The exponential average of W (rather than Wo) follows from similar manipulations: 

(e~ 0W ) - Jdz P^(z ;X )e-^^ 
1 



Z(X ) 

1 

zjxo) 



dz p -l3H(z Tl X T ) _ Z{X T ) 

Z(X ) 



-PAF 



(23) 
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We have used W(zq) = H(z t (zq); X T ) — H(zq;Xo) (Eq, 1 15b)) to get from the first line to the second, and a change of 
variables, z — > z T , to get to the third. 

Eq.l22"lwas originally obtained by Bochkov and Kuzovlev [H, EH, [H], [HI] , whereas Eq. [23] is the nonequilibrium work 
theorem of Refs. 0,0]. These results apply to two physically distinct quantities, Wo and W, corresponding to different 
conventions for defining the internal energy of the system. In each case the exponential average of work reduces to a 
ratio of partition functions. In Eq.[22]the ratio is Z(Xq)/Z(Xq), i.e. unity; while in Eq.[23]it is Z(X t )/Z(Xq), which 
yields the free energy difference AF. 

Let us now consider the special case in which the external forces vanish both at t = and at t = t: 

X Q =X T = (Q,O r --). (24) 

This corresponds to a cyclic process, for which the Hamiltonian begins and ends at H . In this case we have, identically, 
W = Wo fEq. fTBl) and AF = (Eq.[5]). Thus, Eqs. |2"21 and |2"31 are equivalent when the Hamiltonian is varied cyclically. 

Finally, it is instructive to consider a process during which the external forces are switched on suddenly at t = 0, 
from X = (0, 0, • • • ) to X T = (Xi,X 2 , ■ ■ ■ ). Since the process occurs instantaneously (r — > 0), the system has no 
opportunity to evolve, hence z T — zq. Thus, Eq. [TBI gives us 

Wo=0 , W = AH(z ), (25) 

where AH(z) = H(z;X T ) — H(z;Xq), Eq. [22] is immediately satisfied, and Eq. [23] reduces to Zwanzig's perturbation 
formula |20j ] - 

<e-" Aff > = e~ pAF , (26) 
where (• • • )o denotes an average over microstates sampled from the X = (0, 0, • ■ ■ ) canonical distribution. 



IV. EXAMPLE 



Let us now illustrate the general analysis presented above, using the example of a one-dimensional harmonic 
oscillator perturbed by a uniform external force. We take the bare Hamiltonian 

Ho{z) = ^P 2 + (27) 
2m 2 

and we consider a perturbation 

-XQ(z) = -Xq. (28) 

Thus, H = Hq — Xq. The perturbation describes a force X acting along the direction of the coordinate q. The 
canonical distribution at a given force X is 

P cq (z;X) = ^expl-(3(H Q -Xq)}, (29) 

and by direct evaluation of Eq. [S] we get 

F(X) = F(0) - (30) 
Now imagine a process during which the perturbing force is linearly ramped up from zero to some positive value x- 



X t = — , < t < t. (31) 
r 



To simplify the calculations below, we take r to be the period of the unperturbed oscillator: 



r = ^. (32) 
to 



The evolution of the system satisfies Hamilton's equations, 



q = — , p — —mu q -\ , (33) 

m t 
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p(W) = 5(W-AF) 




which can readily be solved. For initial conditions {qo,Po), we get a trajectory 

Qt — <Zo cos(a;i) 4- sm(ujt) H ^-r— [wt ~ sin(a;t)l (34a) 

raw mu 6 T 

X 

Pi = Po cos(o;i) — mujqo sm(ujt) H — [1 — cos(wi)] , (34b) 

hence 

<3V = <7cH o ' Pt=Po- (35) 

mw z 

The quantities Wo and W then follow from Eq. [T5j 

W = x9o-AF , W = AF, (36) 

where 

AF = F( X ) - F(0) = (37) 

From Eq. [36] we obtain explicit expressions for the distributions of work values, po(Wo) and p(W), assuming initial 
conditions (qo,Pa) sampled from equilibrium. Since W — AF for every realization, we have 

p(W) = S (W - AF) . (38) 

In turn, since Wo is a linear function of qo (Eq. I36[) . which is sampled from a thermal, Gaussian distribution with 
mean (q Q ) = and variance cr 2 o = (tow 2 /?) -1 , it follows that Wq is also distributed as a Gaussian, with mean and 
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variance (W ) = -AF, a^ o = X 2 ^ - Explicitly, 

Po(W ) = 



muj 2 f3 



(39) 



It is now straightforward to verify by inspection and Gaussian integration that Eqs.[T]and[2]are satisfied: 

(e-W) = J dW P (W)e-^ w =e-P AF (40) 
(e-? w °) = [ dW oPo (W )e-P w ° = 1. (41) 



The very simple expressions obtained above for Wq and W are consequences of our choice for r, Eq. 1321 The model 
remains solvable for arbitrary r; in that case both p and po are Gaussian distributions, satisfying Eqs. Q] and [H but 
the expressions for their means and variances are more complicated. 




1 
q [arb. units] 



FIG. 2: The bare harmonic potential Uo is shown, along with the distributions of initial particle positions (dark gray) and final 
particle positions (light gray). As we turn on the force X, we shift the distribution by an amount Aq = x/ muj ' 2 — 1-0 away 
from the minimum of the potential, resulting in an average increase in the value of Uo- 



Fig. [T] depicts the distributions po and p. Note that W is negative, while the mean value of Wo is positive. We can 
understand this sign difference as follows. Suppose we adopt the exclusive convention and take the internal energy of 
the system to be given by Hq. We imagine that the particle evolves in a fixed harmonic well 

U Q (q) = to^VA (42) 

under the influence of a time-dependent external force X t . The initial position qg is sampled from a Gaussian 
distribution centered at q = 0, and as we turn on the perturbing force from to x, we displace the particle rightward 
by a net amount Aq = q T — qo = x/ muj2 (Eq. I35|) . The final condition q T is then distributed as a Gaussian whose 
mean no longer coincides with the center of the fixed harmonic well, but rather has shifted by a distance Aq, as 



8 




-0.6 L 



q [arb. units] 

FIG. 3: Same as Fig. [2] only now we think in terms of a time-dependent potential well rather than an externally applied force. 
From t = to t = r the minimum of the harmonic well moves along the path depicted by the thick dashed line. During a given 
realization, the particle suffers a net displacement q T — qo = Aq, hence W = U T (q T ) — Uo{qo) = AF. 



shown in Fig. [2] In effect, the perturbation pushes the particle distribution rightward along the q-axis, and "up" the 
quadratic potential energy landscape, resulting in a positive value for the average work, (Wo) > 0. 

Now suppose that we instead choose the inclusive convention and use the full Hamiltonian H = Hq — Xq to define 
the internal energy of the system. Thus we imagine a particle moving in a time-dependent potential, 

muj 2 ( Aq \ 2 . t 2 , . 

U t (q)=U (q)-X t q= — (q-^-t) + ■ (43) 

where Aq = x/mw 2 and AF — — x 2 /2mu> 2 , as above. Eq. 1431 describes a harmonic well that moves rightward along 
the q-axis with a velocity Aq/r, and slides downward in energy, as depicted in Fig. [3] We can now appreciate 
why W = AF for every realization of the process. From t = to t = r the particle moves by a net amount Aq; 
simultaneously, the well shifts by the same amount, and acquires an energy offset AF: 

Ur{q) = ^f{q-Aq) 2 + AF. (44) 

The particle thus ends with the same displacement relative to the minimum of the well as it began with, so the net 
change in its energy is just the offset A_F. 

In summary, in the exclusive framework, we picture a particle that is pushed rightward by an external force in a 
fixed harmonic well ((Wo) > 0); while in the inclusive framework we imagine a particle that is dragged rightward in 
space and "downward" in energy by a moving harmonic well (W < 0). 



V. WEIGHTED DISTRIBUTIONS 



Here we sketch an alternative derivation of Eqs. Q] and[H 
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Consider an ensemble of realizations of the process described in Section ITT] Let us picture this ensemble as a swarm 
of Hamiltonian trajectories evolving in phase space, represented by a density 2 



which satisfies the Liouville equation, 



f(z,t) = {s[z-zt)), (45) 



= (46) 



Here we use the convenient Poisson bracket notation: {A, B} = (dA/dq) ■ (dB/dp) — (dA/dp) ■ (dB/dq) . In general, 
Eq. |46]does not yield a simple solution; the evolution of f(z, t) can be very complicated, particularly if the underlying 
Hamiltonian dynamics are chaotic. 
For a given trajectory z t , let 



wo(t)= f dt'Y, x k{t')Qk{z t ') (47) 
J o ft 

denote the amount of work performed on th e sy stem to time t, using the definition of work corresponding to the 
exclusive framework (Eq. HOal and Refs . fl2l . fl3l . [Til HH). Since the rate of change of the observable Qk along a 
trajectory z t is given by Qk = {Qk,H} [22], we can rewrite Eq. l47l as 

w {t)= f dt'{X-Q,H} = f dt'{X-Q,H }. (48) 
Jo Jo 

The last equality follows from the identity {X ■ Q, X ■ Q} = . Now consider a weighted phase space density 

g (z,t) = (5(z-zt) exp[-fiwo(f)]j, (49) 
in which each trajectory carries a statistical weight, exp[— (3wo(t)] (see the discussion below). This density satisfies 

?j£ = {H,go}-0{X-Q,Ho}go, (50) 

where the second term on the right accounts for the evolving statistical weights. The derivation of this equation is 
very similar to those found in Section II of Ref. @ and Section 4.1 of Ref. [2lj | . 

Since wq(0) — identically, and since we assume our ensemble is initially prepared in equilibrium, we have go(z, 0) = 
/(z,0) = P eq (z; Xq). Given these initial conditions, the unique solution of Eq. [SDlis the time- independent distribution 

9o(z,t) = — i— exp[-0H o (z)] = P cq (z;X ). (51) 

To see this, note that 

{H, e-^ Ho ] = -P{H, Ho} e-P Ha = /3{X ■ Q, Ho}e~P H °, (52) 

using the derivative rule for Poisson brackets, and the identity {Ho, Ho} = 0. From this result it follows by inspection 
that Eq. ED satisfies Eq. |50j 

The functions f(z,t) and go(z,t) arc two different statistical representations of the same ensemble of realizations. 
Continuing to picture this ensemble as a swarm of trajectories evolving in phase space, / (Eq. |4"S"|) can be viewed 
as a number density, which simply counts how many trajectories are found in the vicinity of z at time t; while go 
(Eq. I49|) can be interpreted as a mass density, if we imagine that each realization carries a fictitious, time-dependent 
mass exp[—(3wo(t)]. Eq. ED then has the following interpretation: when the initial conditions are sampled from 
equilibrium, the "mass density" of trajectories remains constant in time, even as the "number density" evolves in a 
possibly complicated way. Thus while the number of trajectories found near a given point z changes with time, these 



2 See Ref. [2ll for a brief discussion of the ordering of limits implied in Eqs, 1451 and 1491 
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fluctuations are balanced by the evolving statistical weights (fictitious masses) of those trajectories, in precisely such 
a way as to keep the local mass density constant. 

We can obtain analogous results in the inclusive framework (Eq. HObl and Refs. 0, Q). Introducing 

w(t)=- f dt 1 y^X k {t')Q k (z t ,) = - f dt'X-Q, (53) 
Jo k Jo 

along with the corresponding weighted density 

g{z,t) = (&{z-zt) exp[-/M*)])> (54) 

we obtain the equation of motion [H, [2l| 

^ = {H,g}+(3XQg. (55) 

For initial conditions g(z,0) = f(z,0) = P cq (z; Xq), the unique solution is 

9(z, t) = exp[-(3H (z; X t )) = ||g P^(z; X t ). (56) 

The weighted density is no longer constant in time (as was the case with go), but rather is proportional to the 
equilibrium distribution corresponding to the current value of the parameters X. 
The results just obtained are summarized as follows: 

(S(z-zt)exp[-0wo(t)]) = P cq (z;X ) (57a) 

(* (z - z t ) exp[-/Mi)]) = ||g P^(z; X t ). (57b) 

Eqs. [T] and [2] now follow immediately by evaluating Eq. [57] at t = r and integrating both sides over phase space. 
While the derivations presented here arc less elementary than those of Section HTT1 we ultimately gain a stronger set of 
results. By a simple trick of statistical reweighting, we transform an equation of motion that we cannot solve fEq.!46|) 
into one that is easily solved (Eq. [50] or [55]) . The result, Eq. [57] allows us to reconstruct equilibrium distributions P eq 
using trajectories driven away from equilibrium. 

Eqs. I57al and [57bl are in fact equivalent. Multplying both sides of Eq. [57a] by exp[+(3X t ■ Q(z)) and pulling this 
factor inside the angular brackets, we obtain Eq I57bl Conversely, multiplication by exp[— (3X t ■ Q(z)} leads us from 
Eq. I57bl to Eq. I57al However, this equivalence is lost once we integrate over phase space: Eqs. [1] and [2] do not imply 
one another. 

Eq. l57bl can be viewed as a direct consequence of the Feynman-Kac theorem; this observation by Hummer and Szabo 
serves as a starting point for their method of reconstructing equilibrium potentials of mean force from single-molecule 
manipulation experiments carried out away from equilibrium 7]. Moreover, Eq. I57al is essentially a special case of 
Eq. 12 of Ref. 7] (with Wq as generalized by Eq. [551 below), if we assume that their confining potential is initially 
turned off: u(z, 0) = 0. For an alternative approach to estimating potentials of mean force from similar experiments, 
see the "clamp-and- release" method proposed by Adib [23| . 



VI. DISCUSSION 



The nonequilibrium work theorem, Eq. [TJ has generated interest (and controversy [pi 
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26[) primarily for two 



reasons. First, along with the fluctuation theorem for entropy production [271. l28l. 123. l30l. l3ll. |32| . it is one of relatively 
few equalities in statistical physics that apply to systems far from thermal equilibrium. Note that the term "fluctuation 
theorem" has also been used to specify a relation between the response of a system to external perturbations, and a 
correlation function describing fluctuations of the unperturbed system 33, 34]. Second, Eq.[T]predicts that equilibrium 
free energy differences can be determined from irreversible processes, counter to expectations that irreversible work 
values can only place bounds on AF 35] . Eq. [2] shares the first feature - it remains valid far from equilibrium - but 
not the second; it does not seem to be the case that AF can be determined solely from a distribution of values of Wq . 

A crucial distinction in this paper has been the difference between the quantities W and Wo- The recognition that, 
in the literature, various meanings are assigned to the term work, might at first come as an unwelcome surprise. Work 
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is a concept of such central importance in thermodynamics that it ought to be unambiguously defined! However, in 
dealing with a physical situation that involves the mechanical perturbation of a system, the perturbation is inevitably 
accomplished by coupling externally controlled variables (Xi, X%, • ■ • ) to generalized system coordinates (Qi, Q2, ■ • • )■ 
This coupling is represented by a term of the form — ^ fe XkQk (or a nonlinear generalization thereof, see below) in the 
full Hamiltonian that governs the evolution of the system and its surroundings. We are then faced with the question 
of whether or not to view this term as part of the internal energy of the system of interest. As stressed in this paper, 
either choice is perfectly acceptable - this is a question of book-keeping rather than principle - but it is precisely this 
freedom that leads to the ambiguity in the definition of work. For related discussions of this issue, particularly in the 
context of interpretation of experimental data, see Refs. [il l36l. l37l l38l . 

Throughout this paper it has been assumed, following Refs. jl2l. Il3l ILH [TEl]. that the coupling between the forces 
X and the observables Q is linear: H = H — X ■ Q. However, as already observed by Bochkov and Kuzovlev, this 
assumption can easily be relaxed. Had we written the Hamiltonian as 

H(z;X) = H (z)-h(Q;X), (58) 

and assumed h(Q; Xq) = 0, then the entire analysis leading to Eqs. [I] and [2] would have remained valid, provided the 
following definitions of work: 

We recover Eqs. [51 and [TT1 with a linear perturbation h — X ■ Q. 

While the analysis here has been carried out using Hamiltonian dynamics, the conclusions remain valid under other 
frameworks for modeling the evolution of the system. The connection to the stochastic approach taken in Ref. [3] has 
already been noted. Moreover, Eqs. 33 and 34 of Ref. (37j . derived for inertial Langevin dynamics, are equivalent to 
Eqs. [1] and [5] of the present paper. For non-inertial (overdamped) L ang evin dynamics, similar results follow directly 
from the Onsager-Machlup expressions for path-space distributions |39l. |40|. 

Finally, recall the Crooks fluctuation theorem [f| , 

^^=exp[P{W-AF)}, (60) 

where the subscripts refer to two thermodynamic process (forward and reverse) that are related by time-reversal 
of the protocol used to perturb the system. The Bochkov-Kuzovlev papers contain results that are reminiscent of 
Eq. [60l for instance Eq. 7 of Ref. and Eq. 2.12 of Ref. [lij]. However, while Crooks uses a definition of work 
corresponding to W of the present paper, Bochkov and Kuzovlev use Wo, and their results do not involve AF. 
Moreover, in Refs. 

the derivations seem to assume that the initial conditions are sampled from the 
same, unperturbed equilibrium distribution for both the forward and the reverse process (see e.g. Eq. 2.6 of Ref. [3). 
Crooks, by contrast, assumes that the forward and reverse processes are characterized by different initial equilibrium 
states. It would be useful to clarify more precisely the relationship between these sets of results. 
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